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IS.  ABSTRACT 


A  theoretical  investigation  of  thermal  lensing  in  infrared  windows  is 
presented  which  treats  aberration  effects  to  all  orders  in  the  small  angle-of- 
deviation  approximation.  The  model  is  applied  to  a  truncated,  Gaussian, 
infrared  laser  beam  incident  on  a  semitransparent,  isotropic,  disc-shaped 
window.  It  is  shown  that  window  aberrations  limit  the  time  a  diffraction- 
limited  focus  can  be  held  in  the  far-field.  This  diffraction-limited  time  tj  is 
computed  for  some  candidate  window  materials  and  their  relative  merits  are 
discussed.  Some  approaches  to  solving  the  thermal  lensing  problem  from  both 
an  engineering  and  a  materials  point  of  view,  as  well  as  some  program 
research  and  development  needs,  are  discussed. 
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Abstract 

A  theoretical  investigation  of  thermal  lensing  in  infrared  windows  is 
presented  which  treats  aberration  effects  to  all  orders  in  the  small  angle-of- 
deviation  approximation.  The  model  is  applied  to  a  truncated,  G.'.ussian,  infrared 
laser  beam  incident  on  a  semitransparent,  isotropic,  disc-shaped  window.  It  is 
shown  that  window  aberrations  limit  the  time  a  diffraction-limited  focus  can  be 
held  in  the  far-field.  This  diffraction-limited  time  td  is  computed  for  some 
candidate  window  materials  and  their  relative  merits  are  discussed.  Some 
approaches  to  solving  the  thermal  lensing  problem  from  both  an  engineering  and  a 
materials  point  of  view,  as  well  as  some  program  research  and  development 
needs,  are  discussed. 
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Thermal  Lensing  in  Infrared  Laser 
Window  Materials 


1.  INTRODUCTION 

Recent  advances  in  high-power  infrared  lasers  (Emmett,  1971)  suggest  the 
possibility  of  focusing  large  power  densities  in  the  far-field.  If,  however,  a  laser 
beam  of  ronuniform  intensity  passes  through  a  semitrauc ,  arent  window  the  beam 
will  be  distorted  as  the  window  heats  up. 

Distortion  of  a  laser  beam  passing  through  a  liquid  cell  has  been  observed  and 
attributed  (Gordon  et  al,  1965;  Tien  et  al,  1965)  to  a  lensing  action  in  the  liquid 
caused  by  nonuniform  heating.  Theoretical  treatments  of  some  problems  quite 
close  to  the  one  examined  in  this  paper  have  also  been  given  (Tien  et  al,  1965; 
Quelle,  1966;  Foster  et  al,  1970).  In  fact.  Sparks  (1971)  showed  that  a  significant 
amount  of  thermal  lensing  will  occur  in  an  infrared  window  well  before  it  fails  due 
to  other  causes.  For  a  discussion  of  other  modes  of  failure  see  Horrigan  et  al 
(1969).  These  treatments  have  been  restricted  to  a  parabolic  dependence  of  the 
index  of  refraction  with  position,  hence,  no  aberration  of  the  beam  occurs.  Ring¬ 
like  interference  patterns  have  been  observed,  however,  and  attributed  to  aber¬ 
ration  effects  produced  by  thermal  lensing  in  low-loss  liquids  (Whinnery  et  al, 
1967). 

In  this  paper  we  extend  previous  treatments  of  thermal  lensing  (Quelle,  1966; 
Foster  et  al,  1970;  Sparks,  1971)  to  include  window  aberrations  to  all  orders  and 
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we  apply  the  analysis  to  -e  specific  problem  of  calculating  the  amount  of  time  the 
far-field  focused  spot  ru\>ains  diffraction-limited  (Jasperse  etal,  1972).  This 
diffraction-limited  ti.-r.'.-  tj,  turns  out  to  be  independent  of  the  beam  radius,  the 
distance  f.—cc  the  window  to  the  spot,  and  the  Gaussian  focal  lengths  of  the  window. 
Our  analysis  provides  the  laser  systems  engineer  with  one  criterion  for 
choosing  among  *  ‘udidate  windows  and  also  suggests  to  the  materials  engineer 
several  possibllin-.s  for  designing  new  materials. 


2.  A  SOLVABLE  MODEL  FOR  THERMAL  LENSINC  IN  ISOTROPIC. 
DISC-SHAPED  HINDOOS 


2.1  Outline  of  the  Calculation 


K 


- 


In  this  section  we  generalize  previous  treatments  of  thermal  lensing  of  a 
TEMqo  laser  beam  to  include  all  orders  of  aberration  produced  by  a  time- 
dependent  temperature  distribution  in  an  isotropic,  disc-shaped  exit  window. 
Employing  cylindrical  coordinates  (p,$,z)  we  assume  that  the  temperature  dis¬ 
tribution  so  generated  is  independent  of  Furthermore,  we  assume  a  small 
angle  of  deviation  for  a  normally  incident  ray,  allowing  us  to  neglect  the  tempera¬ 
ture's  dependence  on  z  since  it  will  have  a  negligible  effect  on  the  bending  of  the 
ray.  This  means  that,  in  effect,  we  are  dealing  with  a  lens  system  which  has  a 
long  but  finite  focal  length.  The  total  angular  deviation  0W,  of  any  ray  in  such  a 
window,  is  the  sum  of  two  parts: 


0=0  ,  +  0.  , 

w  grad  n  bulge 


(1) 


where  n  is  the  contribu'  '.on  due  to  a  gradient  in  the  refractive  index  n,  and 

©bulge  *s  ttle  contribution  due  to  a  bulging  of  the  originally  parallel  faces  of  the 
window.  Such  an  equation  holds  for  both  the  o-  and  ^-polarized  rays. 

2.2  Calculation  of  n 

We  consider  a  temperature  distribution  in  the  window  having  the  following 
particular  form  when  expressed  in  cylindrical  coordinates: 


T  (p/PD.t)  =  T  (0.0)  +  g(t)T(p/p0)  , 


(2) 


where  T(0, 0)  is  a  constant,  g(t)  is  any  function  of  time,  p  is  the  radial  coordinate, 
and  pQ  is  the  window  radius.  The  term  T (p/pQ)  is  any  monotonically  increa^ng  or 
decreasing  even  function  of  p  expressible  as  the  following  convergent  power  series 

in  p/pq: 
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T (p/pD)  - 


(p/Pq) 


The  T2/  are  dimensionless  constants,  with  TQ  =  1  . 

A  fundamental  theory  accurately  establishing  the  index  n,  in  terms  of  fre¬ 
quency,  temperature,  and  positions  of  all  ions  or  atoms  within  the  solid,  is  not 
available.  In  place  of  such  a  theory,  we  develop  a  phenomenological  treatment  of 
thermal  lensing  by  regarding  n  at  a  given  frequency  as  a  function  of  temperature 
and  the  appropriate  strain  tensor  components.  For  the  case  of  a  plane  stress  in 
a  disc  of  isotropic  material,  the  only  nonvanishing  strain  components  in 
cylindrical  coordinates  are  e  ,  and  eZ2  .  Thus,  in  our  case,  the  refractive 
index  may  be  determined  by  specifying  the  values  of  three  quantities  n.,  i  rep¬ 
resenting  p,  4>,  and  z,  where 


ni  *  ni (T'  £pp*  ur  Uz  > 


For  a  fixed  time,  the  temperature  and  the  strains  depend  only  on  p.  Therefore, 
for  a  small  change  in  n.,  and  following  the  development  given  by  Quelle  (1966)  and 
Sparks  (1971)  we  write: 


In  matrix  form,  these  equations  become: 


3"»/ST\  /  >/%> 
8”z/8T/  y*« 


in  which  np,  n^,  and  nz  are  the  components  of  refractive  index  along  the  three 
orthogonal  directions,  and  2lnp»npnzl*  pp*  ezz)  is  the  "Jacobian-like"  matrix 
of  the  transformation  between  the  n.  and  the  e„,  defined  as: 
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In  an  isotropic  disc  under  plane  stress  and  at  a  fixed  temperature  the  Pockets 
;to-optic  consts 
(Krishnan,  1958): 


elasto-optic  constants  p^,  also  called  the  strain-optic  coefficients,  are  defined  by 


l**\ 


_2 


J. 

.2 


p12 

Pm\ 

1 

P12 

1 

Pll 

1 

P12 

J 

\  P12 

F12 

Pi]/ 

/'pp\ 


8W 


(7) 


nQ  being  the  unstrained  refractive  index.  Under  the  conditions  when  the  n.  are  not 
too  much  different  than  n  ,  each  element  on  the  left-hand  side  of  Eq.  (7)  can  be 

°  O  ___ 

approximated  by -2 (nj-^J/n^  .  This  allows  Eq.  (7)  to  be  rewritten  as: 


n  -  n  \ 
P  0 

3 

fpll  p12  p12 

f'pp\ 

-n 

ntf'no 

0 

■  2 

P12  P11  p12 

nz"noj 

1P12  Pl2  Pllj 

iczz  j 

(8) 


With  the  aid  of  Eqs.  (8),  one  can  calculate  each  9n./9e..  so  that  the  J -matrix  in 
Eq.  (6)  can  be  replaced  by: 
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The  thermally  induced  stress-strain  relations  are  (Boley  et  al,  1960): 


(10) 


where  Y  is  Young's  modulus,  the  ...  are  the  three  stress  components  in  cylindrical 
coordinates,  v  is  Poisson's  ratio,  and  a  is  the  linear  expansion  coefficient.  Nov/ 
we  can  define  the  initially  unstrained  "fiducial  state"  as  that  for  which  all  o„  are 
zero  at  zero  time  and  at  the  uniform  temperature  T(0, 0).  Then,  any  change  in 
temperature  and/or  stresses  produces  the  following  change  in  strains: 


(11) 


From  Eqs.  (6)  note  that  each  dnjBT  is  measured  at  constant  strain.  Since  each 
n.  is  not  too  much  different  from  nQ,  there  would  be  no  loss  in  generality  if  each 
3nj/9T  were  replaced  by  the  derivative  evaluated  in  the  fiducial  condition,  namely, 
(3n/3T]n~n  .  Substituting  this  information,  together  with  Eqs.  (9)  and  (11),  into 
(6),  and,  using  the  fact  that  the  &n.  can  also  be  expressed  as  (nj-nQ),  we  get  after 
simplifying: 


where 


with  n  denoting  the  average  index  of  refraction, 

CoS  Pll*2,/Pl2  * 


Cl2  ^Pw  '  "Pll  * 

,  -  •* 
and 


(13) 


(14) 

(15) 


AT  =  T  (p/po.t)  -  T  (0,0)  . 


(16) 


The  C's  are  related  to  the  parallel  and  normal  components  of  the  stress-optic 
coefficients  (B^  and  B^,  respectively)  via: 


B 


II 


(17) 


and 


B. 


nl  c 

2Y  U1  * 


(18) 


For  plane  stresses  in  isotropic  materials  with  disc-type  geometry,  the 
relationship  between  the  stress  and  the  temperature  change  from  the  fiducial  state 
is  (Boley  et  al,  1960): 


MMiilWliBiiia 
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The  equation  for  n  is  not  written  out  since  its  effect  on  the  bending  of  a  ray 
z 

in  the  small  angle-of-deviation  approximation  is  negligible. 

The  fact  that  the  indices  of  refraction  are  different  for  each  polarization  shows 
that  the  window  develops  a  stress-induced  birefringence. 

2.3  Calculation  of  dgrad  n 

Consider  a  ray  normally  incident  at  a  distance 
p.  from  the  z  axis  of  a  parallel-faced  disc,  having  an 
index  n(p/p0,t)  as  shown  in  Figure  1.  A  nonlinear, 
second-order,  differential  equation  for  the  ray's 
trajectory  (p  as  a  function  of  z)  in  the  small  angle-of- 
deviation  approximation  is  solved  approximately  in 
Appendix  B.  Its  solution  is: 

psp.  [l+Kz2  +  oj|K|2z4j],  (22) 

where 

Apply Snell's  Law  at  the  exit  surface  of  the  window, 
we  see  that: 


z=L„ 
o 

where  Lq  is  the  initial  thickness  of  the  window.  From  Eq.  (22),  including  up  to 
second  order: 

7i£  =2PiKL0.  (25) 

z=Lo 

Utilizing  Eqs.  (23)  and  (25)  in  (24),  we  obtain  a  stress-induced  birefringence  in 
®grad  n  given  for  p  and  4>  polarizations  by: 


n(0  t)  =*  .greeted  anJ 
~  incident  angl 


\ 


R« 


Figure  1.  The  Bending 
of  a  Ray  as  It  Passes 
Through  a  Disc-Shaped 
Window 
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(26a) 


grad 


X  2L  g(t)  oo  ,  /  p,  \2I-1 

n.2PlLo„(0.«K *.-f—  . 


(26b) 


Throughout  this  paper  we  adopt  the  convention  that  the  angle  that  any  deviated 
ray  makes  with  the  horizontal  will  be  positive  if  the  ray  bends  away  from  the 
z  axis  and  negative  if  towards  the  z  axis. 


2.4  Calculation  of  flbulge 


The  equation  for  the  plane  curve  L  associated  with  one  face  of  the  bulging 
window  is: 


L(p/p0.t)  =  Ln+TLrt  A  e„  . 


I  '  n  U  C__ 

o  z  o  zz 

Utilizing  Eqs.  (11)  for  Ac_  .  Eqs.  (A-4)  and  (A-5)  for  o  and  cr , ,,  and 


(27) 


'zz'  -'1"  - - * - -  —  ’pp  -  $$ 

Eq.  (A-l)  for  AT,  the  above  equation  can  be  rewritten  as: 


L  g(t)  a  (1+  v)  oo 

L(p/po.t)*L(0.t)+  g -  S  T?|  (p/p0)  21 

I  “1 


(28) 


in  which,  as  usual,  all  the  terms  independent  of  p  have  been  collected  under  the 
symbol  L  (0,  t).  In  Appendix  C,  we  show  that  6Kl|1p0  can  be  expressed  as: 


0,  .  =2 

bulge 


fa"’(£L 


(29) 


'P=Pi 

Utilizing  Eq.  (28)  results  in: 


0, 


/2L  g(t)\  oo  t2f-l 

w  *(-£-7*  t2<  1  w 


(30) 
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2.5  The  Window  Angle  of  Deviation,  &w 

The  total  angle  of  deviation  due  to  the  window  itself,  0^,  (where  j  can  stand 
for  either  p  or  ^  polarizations)  can  be  obtained  by  inserting  Eqs.  (26  a  or  b)  and 
(30)  into  (1).  It  is  * 


+  a  (l+v)(n-l)  T2i 


m 


2/  —  1 


(31) 


The  Gaussian  focal  lengths  R^  are  related  to  the  terms  of  0^,  which  are  linear 
in  p..  They  are 


•D  3 

Rg-8J 

w 


2  LQg(t)|  n32  +  « (l+t/)(n-l)  T2f 


(32) 


We  note  that  all  formulas  in  Section  2  apply  for  both  transient  and  steady-state 
conditions,  since  g(t)  may  be  any  function  of  time. 


3.  APPLICATION  TO  A  TRUNCATED  GAUSSIAN  LASER  BEAM  INCIDENT 
ON  A  DISC-SHAPED,  ISOTROPIC  UNDO* 


3.1  Assumptions 

In  applying  the  solvable  model  given  in  Section  2  to  the  specific  problem  of  a 
laser  pulse*”'  of  time  t  incident  on  a  semitransparent,  disc-shaped  window,  we 
make  the  following  assumptions: 

-o 

(1)  total  angular  deviations  are  less  than  about  10  radians; 

(2)  the  window  material  is  isotropic; 

(3)  the  thermally  induced  stresses  and  the  corresponding  strains  are  per¬ 
fectly  elastic; 


=3 

M 

13 

-  3 

il 

t; 

I 

h 


*  The  first  two  terms  of  each  of  these  series  have  also  been  derived  by  a 
Raytheon  group  using  a  different  method,  namely,  by  computing  phase  differences 
between  rays  traversing  the  window  at  different  radial  distances.  Their  phase 

angles  are  related  to  our  0^,  via  f  k  0^  dp.  See  F.  A.  Horrigan  and  T.  F.  Deutsch, 

w  w 

o 

Raytheon  Research  Division,  Second  Quarterly  Report,  Contract  No.  DAAH01-70- 
C-1251  (1970). 

**  The  mode  of  laser  operation  contemplated  here  is  interruoted  cw  with 
typical  values  of  t  on  the  order  of  a  few  seconds. 
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(4)  the  laser  pulse  time  is  less  than  or  comparable  to  the  smallest  of  the 
characteristic  times  for  heat  flow; 

(5)  bulk  absorption  dominates  other  loss  mechanisms  in  the  calculation  of  the 
radial  temperature  gradient;  and 

(6)  the  power  per  unit  area  absorbed  in  the  window  PA  is  given  by 

PA(p/Pq)  =  PI  (P/Po)  (1'R)  (1_e  ^AL°)(l-Re  PaL°  )  •  (33) 

Here  Pj  is  the  incident  power  per  unit  area,  R  is  the  reflectivity  of  the  window 
material,  and  £A  is  the  bulk  absorption  coefficient. 

For  0aLq  «  1,  Eq.  (33)  reduces  to: 

PaWpo^WpoMSaV  (34) 


3. 2  Calculation  of  the  Temperature  Distribution 

The  equation  of  heat  flow  in  the  window  is: 


(35) 


where  K  is  the  thermal  conductivity,  c’  is  the  specific  heat  times  the  density,  and 
Q,  given  by  Pa/Lq,  is  the  rate  of  heat  generated  within  the  window  per  unit  vol¬ 
ume.  To  give  the  reader  some  appreciation  of  the  temperature  distribution  in 
cylindrical  geometry,  two  model  problems  are  presented  from  Carslaw  and 
Jaeger  (1959). 

(1)  Consider  a  disc,  infinite  along  the  p  dimension  and  finite  in 

z  (-Lq/2SzSL0/2)  .  With  zero  initial  temperature,  zero  surface  temperatures, 

and  Q  constant  for  t  >  0,  the  solution  for  T(z,  t)  is  given  by  Eq.  (7),  page  130  of 

Carslaw  et  al  (1959).  Its  graph,  depicted  in  our  Figure  2,  plots  normalized 

temperature  versus  normalized  distance  with  normalized  time  as  parameters. 

The  dominant  term  in  the  infinite  series,  namely,  the  first  term,  is  proportional 

to  { l-exp(-t/r  )]  ,  where  r  ,  the  characteristic  time  for  heat  flow  along  the 

z  axis,  is  given  by  c*  Lr  /Kv  .  Consequently,  for  t  <  r  ,  we  see  that  to  a  fair 
o  z 

approximation  the  temperature  varies  linearly  with  time. 

(2)  Consider  the  same  disc  as  in  (1)  above.  Besides  a  zero  initial  tempera¬ 
ture  and  constant  Q  for  t  >  0,  both  faces  at  z  =  ±Lq/2  transfer  heat  into  a  medium 


Figure  2.  Normalized  Temperature  vs 
Normalized  Axial  Distance  for  the  Case 
of  a  Disc,  with  Infinite  Radius,  at 
Zero  Initial  Temperature,  Zero  Surface 
Temperatures,  and  Constant  Rate  of 
Heat  Generation.  The  numbers  on  the 
curves  refer  to  normalized  times 
(8t/ir  2T  ) 


0.0  02 


0.4  0.6 
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at  zero  temperature.  The  solution  for  T(z,t,h)  is  given  by  Eq.  (12),  page  132  of 
Carslaw  et  al  (1959).  Its  graphs,  labelled  in  exactly  the  same  fashion  as 
Figure  2,  are  depicted  in  Figure  3  for  three  different  dimensionless  radiation 
parameters,  LQh/2;  h  is  defined  as  the  ratio  of  the  surface  heat  transfer  co¬ 
efficient,  or  surface  conductance,  to  the  conductivity.  Since  the  CGS  dimensions 
of  the  surface  transfer  coefficient  are  cal/cm2-°C-sec,  then  h  has  dimensions  of 
reciprocal  length.  Problem  (1)  could  also  be  considered  as  a  special  case  of 
problem  (2),  in  which  LQh/2  is  infinite. 

Returning  to  our  finite  window,  we  choose  to  represent  the  power  per  unit 
area  of  the  incident  beam  as  any  monotonically  decreasing  even  function  of  p  ex¬ 
pandable  in  a  convergent  power  series  of  the  form: 


Pr  <p/p0>  =  P(0)^S  (“D*  P2i  (p/p0) 


Assuming  time  independence,  and  utilizing  Eq.  (34),  Q  can  be  expressed  as: 
P a  (p/p  ) 

Q<p/p0>  ■  — t  °  “  0APr  <p/p0>  • 


If.  in  the  heat  equation,  the  Laplacian  term  becomes  negligible  compared  with 
Q/c',  then  direct  integration  results  in  T  being  linearly  dependent  on  t.  That  is. 


T(p,  t)  st  Q/c'  +  T(0,0)  . 


(38) 
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Figure  3.  Normalized  Temperature  vs  Normalized  Axial  Distance  for  the  Case  of 
a  Disc  with  Infinite  Radius,  Zero  Initial  Temperature,  Constant  Rate  of  Heat 
Generation,  and  Both  Faces  Transferring  Heat  into  a  Zero  Temperature  Medium 
with  Normalized  Radiation  Parameters  (a)  Loh/2  =  10;  (b)  Loh/2  =  1. 0;  and 
(c)  LQh/2  =  0. 1.  The  numbers  on  the  curves  refer  to  normalized  times  (8t/ir2rz) 


Initial  uniform  temperature  has  been  assumed,  and  the  z  dependence  las  been 
neglected  because  it  has  negligible  effect  on  the  bending  of  a  ray  in  the  small  angle- 
of-deviation  approximation.  Substituting  Eqs.  (37)  and  (36)  into  (38)  produces: 


Pq  must  also  equal  unity. 

K  2  I 

Let  us  now  consider  under  what  conditions  the  Laplacian  term,  V  T(p,  t)  , 

would  be  much  less  than  |  Q/c'  |  .  First,  we  assume  that  T  has  the  form  depicted 
in  Eq.  (39).  T  len.  for  cylindrical  coordinates: 


From  Eqs.  (37)  and  (36): 


^  oo  2  £ 

P(0)  4Jo  (-D1  P2l(%)  •  (43) 

For  those  cases  in  which  both  sums  converge  over  the  range  of  p  and  have  values 
comparable  in  magnitude,  it  follows  that  |  (K/c')  V3t|  «  |Q/c'  |  when  t«c'p3/4K. 

We  designate  this  quantity,  c'p3/4K,  as  t^,  the  characteristic  time  for  heat  flow 
along  the  p  direction.  It  should  be  noted  that,  in  general,  any  characteristic  time 
for  heat  flow  depends  on  cooling  details. 

Thus,  we  see  that  temperature  will  be  approximately  a  linear  function  of  time 

o 

so  long  as  t  <  and  r  ,  respectively.  Since  the  ratio  Tp/Tz  goes  as  (pQ/ L0)  , 
which  is  much  greater  than  unity,  then  t  is  by  far  the  smaller  of  the  two  quantities. 
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(d)  LINEAR  APFROXIMATION- 
Eq  (39) 


-(b)  Loir ,  | 


■OJ^.,0 


as 


Figure  4.  Normalized  Temperature  vs 
Normalized  Time  at  the  Center  of  the 
Infinite  Disc  Referred  to  in  Figure  3. 
Curves  (a)  through  (c)  apply  to  Figures 
3(a)  through  3(c),  respectively,  while 
curve  (d)  is  the  approximate  linear 
solution  given  by  Eq.  (39) 


Consequently,  we  can  conclude  that  the 

temperature  behavior  given  by  Eq.  (39) 

is  approximately  valid  for  times  less 

than  t  . 
z 

Numerical  values  of  tz  for  some 
materials  under  consideration  will  be 
quoted  later. 

We  can  get  a  quantitative  estimate 
of  just  how  valid  this  linear  approxima¬ 
tion  is  by  reconsidering  problem  (2) 
mentioned  above.  Going  back  to  Figures 
3(a)  through  3(c)  we  can  obtain  informa¬ 
tion  about  how  normalized  temperature 
varies  with  normalized  time  for  the 
three  different  radiation  parameters  at, 
say,  the  center  plane  of  the  disc 
(2z/Lq  =  0).  We  plot  this  data  from  these 
three  figures  in  Figure  4  as  curves  (a) 
through  (c),  respectively.  The  line  (d) 
represents  the  approximate  linear  solu¬ 


tion  given  by  Eq.  (39).  From  Figure  4  we  see  that  as  8t/ir  tz  varies  from  0  to  1 
the  approximate  linear  solution  is  within  about  3%  of  the  exact  solution  for 
LQh/2  =  0. 1,  within  about  9%  of  the  exact  solution  for  LQh/2  =  1.  0,  and  within  24% 
for  LQh/2  =  10. 

In  the  case  of  a  finite  disc  in  which  po/LQ  is  still  much  larger  than  unity  and 
for  a  given  h,  we  can  arrive  at  a  more  accurate  g(t)  function  than  the  approximate 
linear  solution.  From  the  appropriate  normalized  temperature  vs  distance  curves, 
as  shown  in  Figure  3,  for  example,  we  could  average  each  time-curve  over  the  z 
coordinate,  thereby  eliminating  the  z-dependence  of  temperature.  From  this  in¬ 
formation,  a  normalized  temperature  vs  time  curve,  similar  to  those  shown  on 
Figure  4,  could  be  plotted  and  a  g(t)  function  could  be  curve- fitted  from  it  for  use 
in  Eq.  (2). 


3.3  Angular  Deviations  for  Truncated  Gaussian  Beams  in  the 
Quartic  Polynomial  Approximation 

If  we  choose  the  incident  beam  as  having  a  truncated  Gaussian  power  distribu¬ 
tion  given  by 


PiWPq)  *  P(0)  exp  [-w  <p/p0)2] 


(44) 
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for  0  £  p  £  pQ,  where  W  is  a  dimensionless  constant,  and  make  a  qua r tic  poly¬ 
nomial  approximation  to  it,  we  obtain 


lyp/p,,)  “  p<o>  [i  -  p2  (p/p0>2  +  p4  <p/p0»4]  - 


(45) 


Performing  a  least  squares  fit  to  determine  P2  and  P^,  we  find  that  for  W  =  1, 

P2  =  0. 96  and  P^  =  0.  33,  and  that  Eq.  (45)  fits  Eq.  (44)  to  within  1%  everywhere; 
whereas  for  W  =  2,  P2  =  1.75  and  P4  »•  0.  93,  and  Eq.  (45)  fits  Eq.  (44)  to  within 
32%  everywhere. 

The  formulas  for  6^  and  0^  in  the  quartic  polynomial  approximation  (that  is, 
1=1  and  2  only)  are  given  by  Eq.  (31)  in  conjunction  with  Eqs.  (21),  (40),  and 
(41). 

They  are: 


((l-3t/)p11+(3-5%)12J  +  a(l+v)(n-l)j 


3  3 

+  2  P4pfl.o  *  T2  1 1+  +  *  <1W><"-  »](^-) 


(46a) 


and. 


w 


'  °  'I  0=0  J  ° 

+  2 p4 +T^((5-v)p11+(l-llv)p12)  +  '«(l+r)(n-l)j^L^  J.  (46b) 


Utilizing  the  substitutions: 


FS<«  P(0)P2t  [(lr)o,0+ T  ( (l-MPn+H-SvIPu)  . 

HI- 


+  a(l+v)(n- 


(17a) 


We  note  that  the  Fj  and  Fg  are  functions  of  material  parameters;  Fj  may  be 
of  either  sign,  but,  for  all  materials  studied  so  far,  Fg  is  always  positive  and  is 
typically  less  than  unity.  Under  these  circumstances,  the  qualitative  behavior  of 
will  be  governed  by  the  sign  of  F^  .  That  '■»,  depending  on  whether  F^  is 
positive  or  negative,  the  window  may  act  like  a  converging  or  diverging  lens, 
respectively,  having  spherical  aberration  for  each  polarization.  If  more  terms 
were  kept  in  the  approximation  to  the  beam  intensity,  higher-order  aberrations 
would  appear. 


4.  THE  DIFFRACTION-LIMITED  TIME, 

Diffraction  patterns  in  the  scalar  (Campbell  et  al,  1969)  or  vector  (Bendow  etal, 
1972)  Kirchhoff  approximations  can  easily  be  calculated  using  modern  computers. 
Such  calculations  do  not,  however,  give  the  materials  engineer  much  insight  into 
how  to  design  a  better  window,  nor  do  they  give  the  laser  systems  engineer  a 
simple  way  of  assessing  the  performance  of  known  window  materials.  The  object 


^SpB3=S5i 


18 


of  this  work  is  to  obtain  approximate  formulas  of  use  to  both  types  of  investigators. 
To  that  end  we  give  an  approximate  treatment  of  diffraction  effects  using  a  combina¬ 
tion  of  geometrical  optics  and  scalar  Kirchhoff  theory  applied  separately  to  the 
p  and  <j>  polarizations. 

Initially,  the  spot  size  at  focus  in  the  far-field  is  determined  only  by  diffrac¬ 
tion.  Aberrations  produced  by  the  window  eventually  dominate  diffraction  effects, 
however,  and  rapidly  enlarge  the  spot  size.  We  define  the  diffraction-limited  time 
td  as  the  time  a  laser-window  system  takes  to  pass  from  its  diffraction-limited 
condition  at  t  =  0  to  its  window -aberration-limited  ■  jndition  at  t  =  t^.  In  this 
section  we  compute  t^  and  show  that  it  depends  on  window  aberrations. 

The  diffraction  half-angle  0^.^  is  a  convenient  approximate  measure  of  the 
amount  of  spreading  of  a  beam  due  to  diffraction  effects.  We  take  it  as  that  half¬ 
angle  which  contains  90%  of  the  total  diffracted  power  at  the  focused  far-field  spot. 
This  can  be  obtained  approximately  by  resorting  to  Figure  4  of  Olaofe  (1970).  A 
horizontal  line  drawn  from  the  ordinate  0.  9  (which  represents  fractional  power) 
intersects  the  q  =  1  curve  (which  corresponds  to  our  W  =  1  case)  at  an  abscissa 
value  of  approximately  3.  This  abscissa  is  the  normalized  radial  distance, 

( 2v/X)pQ0 where  X  is  the  wavelength  of  the  radiation.  Equating  the  latter 
expression  to  3,  we  obtain: 


e 


diff  2-irp 


(50) 


In  the  case  of  a  beam  prefocused  to  have  a  focal  length  Ro,  the  incident  angle 

of  deviation  is  given  by  -p/R  .  It  can  be  added  to  0^  [see  Eq.  (49)]  to  obtain  the 

o  w 

total  angle  of  deviation: 


'!«£[>- ■*(£-)’] 


(51) 


in  which  p  has  replaced  p^  . 

Because  of  the  presence  of  the  aberrating  terms  in  Eq.  (51),  the  rays  incident 
at  different  p  will  not  all  come  to  a  focus  at  the  same  range.  Instead,  the  envelope 
of  these  rays  will  form  a  circle  of  least  confusion.  In  Appendix  D  we  calculate  the 
time  dependence  of  the  radius  and  location  of  this  circle  of  least  confusion  using 
geometrical  optics,  as  depicted  in  Figure  5  for  the  quartic  polynomial  approxima¬ 
tion. 


am 


Figure  5.  Two  Rays  Emanating  from  Two  Different 
Radii  Within  the  Window  for  the  Case  When  7  Increases 
Monotonically  with  p.  Their  intersection  (D,  r)  will 
give  the  distance  and  radius,  respectively,  of  the 
circle  of  least  confusion  after  r  is  maximized 


The  diffraction-limited  time  occurs  when  the  radius  of  the  circle  of  least  con 
fusion  has  grown  so  that  it  is  equal  to  the  radius  of  the  diffraction  circle.  This 
latter  radius,  subtended  by  the  angle  0diff  of  Eq.  (50),  can  be  calculated  in  the 
small-angle  approximation  by  multiplying  fldiff  by  its  distance  from  the  window. 
At  the  diffraction-limited  time  td  this  distance  is  also  equal  to  (td),  as  given 
in  Eq.  (D-ll)  of  Appendix  D.  Thus, 


rc(td)  =  0diff  Dc  V* 


(52 


Substituting  Eqs.  (50),  (D-10),  and  (D-ll)  into  Eq.  (52)  results  in: 


>kl(td>  ?  wr- 


(53 


Depending  on  the  polarization,  either  Eqs.  (47a,  b)  or  (48a,  b)  can  be  employed  in 
the  above  equation.  Solving  for  the  time  gives: 


•s- 


(rir)(pfe) 


c'/0« 


(It)  +%  (a-5 V)p ! ^(5-7^12) +Cf(l+v)(n-1) 


(54a 


for  the  p  polarization,  and. 


(^) 

/  1  \ 

'■/<sA 

\P<0)pJ 

j(ft) 

3 

+22- 

12 

J=0 

'(5-v)pn+(l-llv)p12)+a(l+i/)(n-l) 

(54b) 


for  the  <j>  polarization.  Note  that  tjjj  and  t^  are  independent  of  Ro,  pQ,  the  quadratic 
term  P^,  and  the  Gaussian  focal  lengths  of  the  window  ( see  Eqs.  (32)  and  (21)). 

The  diffraction-limited  time  t^  serves  as  an  excellent  figure  of  merit  for  thermal 
lensing  since  it  is  independent  of  these  quantities.  Note  also  that  the  difference 
between  t^  and  t^  for  a  given  material  serves  as  a  measure  of  the  importance  of 
birefringence  in  determining  the  true  diffraction  patterns. 

Figure  6  gives  a  graphical  interpretation  of  t^  for  the  cases  of  Dq  either  de¬ 
creasing  or  increasing  with  time. 


Figure  6.  The  Diffraction  Half- 
Angle  Sdiff  and  the  Locus  of  the 
Radius  of  the  Circle  of  Least 
Confusion  r  (t)  for  Two  Different 
Cases 


Dc(t)  DECREASING  WITH  t 


There  are  two  other  quantities  of  interest  which  are  calculated.  The  first  is 
the  temperature  difference  between  the  center  of  the  window  and  its  edge,  denoted 
by  AT(t).  Referring  back  to  Eq.  (39)  and  employing  only  up  to  quartic  terms  in 
p/pQ,  (that  is,  1=2)  the  temperature  difference  becomes: 


AT(t)  s  T(0,  t)  -  T(l.t) 


P(0)  (P2-P4)t. 


(55) 


The  second  quantity  of  interest  is  the  amount  of  bulging  in  one  face  of  the 
window,  denoted  by  AL  (t).  Referring  back  to  Eqs.  (28),  (40),  and  (41),  and, 
proceeding  as  in  the  calculation  of  AT  (t)  above,  we  get: 

AL(t)  s  L(0,  t)  -  L(l,t)  =(^f)p(0)(P2-P4)  ^-a(l+v)t  .  (56) 

Thus,  both  AT(t)  and  AL(t)  vary  linearly  in  time  in  our  approximation. 

Values  of  t£,  t^,  characteristic  time  r  ,  as  well  as  those  of  AT  and  AL 
evaluated  at  the  time  t£j  ,  have  been  computed  for  a  number  of  candidate  window 
materials  for  which  reasonably  complete  data  have  been  found.  They  are  shown 
in  Tables  1  and  2.  In  these  computations  we  considered  a  Gaussian  incident  beam 
truncated  at  the  e  power  point  and  chose  Lq  =  3  cm  with  P(0)  =  500  watts/cm  . 
Note  that  for  times  approaching  td  in  the  I  -  VII  and  II  -  VII  compounds  the  window 
experiences  radial  temperature  differentials  from  center  to  edge  in  the  tens  of 
degrees,  while  the  center  of  each  face  bulges  out  by  factors  several  times  that  of 
the  wavelength.  The  corresponding  values  for  the  other  substances  are  very 
much  smaller.  * 

The  material  dependent  parameters  used  in  these  calculations  were  taken 
from  the  literature  and  are  listed  in  Tables  3  and  4.  These  parameters  should  be 
considered  as  being  representative  only  and  are  by  no  means  definitive.  The 
reason  for  this  is  that  oftentimes  the  literature  will  report  either  many  discrepant 
values,  questionable  values  (for  example,  no  listing  of  the  wavelength  or  tempera¬ 
ture  at  which  the  data  was  measured),  or  no  values  at  all.  Data  on  thermal  con¬ 
ductivities  seems  to  be  especially  erroneous.  In  such  instances,  one  has  to  guess 
which  values  seem  more  reasonable  compared  with  those  of  similar  materials, 
or,  to  calculate  an  average  value  either  from  among  all  the  listed  data  or  from 
that  which  seems  reasonable,  or,  to  extrapolate  from  known  data. 

The  value  of  Poisson's  ratio  hardly  changes  from  substance  to  substance  and 
was  taken  to  be  0.  3  throughout.  The  only  exception  to  this  was  found  to  be  TI 

No.  20  glass  for  which  v  is  0. 19. 

-2 

For  an  e  truncated  Gaussian  beam  the  relative  ordering  of  the  materials 
by  t^  and  td  will  not  change  but  all  of  their  values  will  decrease  by  a  factor  of 
0.  33/0. 93  or  by  about  1/3. 

We  emphasize  here  that  all  formulas  in  this  section  are  valid  only  for  times 
less  than  or  comparable  to  tz  . 

*  We  note  here  that  the  diffraction  patterns  for  windows  used  under  actual 
conditions  will  be  very  sensitive  to  small  changes  of  shape  of  the  window  produced 
by  nonthermal  effects,  that  is,  those  produced  by  small  vibrations  or  other 
mechanical  disturbances. 
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Table  2.  Calculations  for  Some  Window  Materials  for  3/i  <,  X  £  5*< 


Table  4.  Values  of  Parameters  Used  in  Calculations  for  3m  £  *  £  5m 


*  Estimated  values 
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S.  DISCUSSION 

5.1  Di(Tr»ctio»-Li«ited  Ti»e 

Equations  (54  a,  b)  show  that  the  material-dependent  portion  of  t^  is  directly 
proportional  to  the  ratio  c'/0A  and  inversely  proportional  to  the  sum  of  three 
terms.  The  first  two  terms  of  this  sum  are  contributions  to  the  spherical  aber¬ 
ration  coming  from  the  gradient  in  the  refractive  index.  We  will  refer  to  the 
second  quantity  as  the  "elasto-optical"  term.  The  third  term  arises  from  the 
window's  bulging,  and,  under  the  prevailing  circumstances,  always  makes  a  posi¬ 
tive  contribution.  With  regard  to  the  substances  studied,  the  numerical  values 
corresponding  to  each  of  these  terms,  including  the  elasto-optic  for  both  polariza¬ 
tions,  are  listed  in  the  last  four  columns  of  Tables  1  and  2. 

For  the  I  -  VII  and  II  -  VII  compounds,  (3n/3T)a=o  is  negative  and  of  com¬ 
parable  magnitude  to  <*(1+ v)(n-l).  Therefore,  these  two  terms  almost  cancel. 

This  allows  the  elasto-optic  terms  to  play  a  very  significant  role  even  though  they 
are  typically  an  order  of  magnitude  smaller. 

For  covalent  crystals  (3n/3T)ff  is  observed  to  be  positive  and  larger  than 
o(l+v)(n-l),  which,  in  turn,  is  larger  than  either  elasto-optic  quantity.  Thus, 
the  elasto-optic  terms  will  have  a  small  effect. 

In  glasses,  (3n/9T)a=o  and  a(l+v)(n-l)  are  both  positive,  comparable  in  mag¬ 
nitude,  and  only  slightly  greater  than  the  elasto-optic  quantities.  The  only  glaring 
exception  is  As2S3  glass,  for  which  the  elasto-optic  terms  at  4p  are  comparable 
to  ar(l+v)(n-l)  and  both  elasto-optic  terms  are  at  least  an  order  of  magnitude 
greater  than  (3n/3T)0=o  . 

One  should  recall  that  the  values  of  t^  given  in  Tables  1  and  2  are  accurate 
only  if  the  assumptions  listed  in  Section  3  are  valid.  Actual  materials  may  not 
perform  as  well. 

3.2  Relative  Merits  of  Candidate  Window  Materials 

One  must  also  remember  that  the  diffraction-limited  time  is  but  one  factor 
in  the  selection  of  a  window  material  for  a  particular  application.  In  the  following 
discussion  of  the  relative  performance  of  materials  listed  in  Tables  1  and  2,  we 
will  include  some  of  these  other  factors. 

For  10.  6-micron  operation  the  alkali  halide  crystals  listed  in  Table  1  are 
much  better  in  their  ability  to  hold  a  diffraction-limited  focus  than  any  other  class 
of  materials  studied.  Their  single-crystalline  mechanical  properties  are,  how¬ 
ever,  generally  poor.  They  are  hygroscopic  and  they  cannot  easily  be  coated  for 
protective  or  antireflective  purposes.  The  II-VI  and  III-V  crystals  have  fair 
mechanical  properties  but  are,  in  general,  only  marginal  in  terms  of  beam 
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distortion.  In  addition,  one  expects  some  difficulty  in  growing  large-size  optical- 
quality  crystals  of  these  materials.  From  a  thermal  lensing  point  of  view  CdTe 
appears  to  be  quite  a  bit  better  than  ZnSe  or  the  III-V's  examined  in  Table  1.  Per¬ 
haps,  as  better  samples  of  the  other  materials  are  prepared,  thereby  producing 
lower  values  of  0A,  for  example,  these  samples  may  become  more  competitive 
with  CdTe.  Germanium  has  a  rather  low  diffraction-limited  time  but  its  other 
properties  are  good  except  for  its  large  index  of  refraction.  With  an  index  of  4, 

Ge  would  probably  require  antireflective  coatings.  The  TI  glasses  have  good 
mechanical  properties  and  can  be  cast  in  large  sizes  but  have  rather  low  values 
for  t^.  In  addition,  they  have  characteristic  thermal  times  (rz  and  r  )  which  are 
quite  large  and  would  probably  require  rapid  face  cooling  even  for  pulsed  laser 
operation.  The  polycrystalline  materials  shown  in  Table  1  have  values  of  td  that 
are  too  low  but  have  good  mechanical  and  thermal  properties.  In  general,  poly¬ 
crystalline  materials  would  show  more  promise  if  techniques  could  be  found  to 
reduce  their  bulk  absorption  coefficients  closer  to  intrinsic,  s’ngle-crystalline 
values. 

The  situation  for  3-  to  5-micron  laser  operation  is  somewhat  better.  The 
I-VH  and  II-VII  compounds  shown  in  Table  2  all  appear  to  be  good  candidate 
windows:  they  hold  a  diffraction  limited  focus  for  a  long  time:  they  have  good 
mechanical  and  thermal  properties,  are  not  hygroscopic,  and,  from  all  indications, 
large  optical  quality  crystals  can  be  grown.  The  remaining  materials  shown  in 
Table  2  are,  at  best,  only  marginal  from  a  thermal  lensing  point  of  view.  As 
more  data  becomes  available,  other  materials  may  also  be  found  which  are  good 
candidates  for  3-  to  5-micron  operation. 

6.  SOME  APPROACHES  TO  SOLVING  THE  THERMAL  LENSING  PROBLEM 

6.1  Muterial-Type  Solutions 

6.1.1  MINIMIZE  0A 

One  of  the  most  obvious  ways  to  minimize  aberration  effects  is  to  find  a 
material  with  a  very  small  0A  and  with  negligible  beam  energy  losses  by  other 
mechanisms.  This  may  involve  synthesizing  a  very  pure  solid  or  designing  a 
completely  new  material. . 

6. 1. 2  MUTUAL  CANCELLATION  OF  ABERRATING  TERMS 

From  an  inspection  of  Eq.  (54)  it  is  clear  that  t^  is  large  if  the  denominator 
in  the  material  dependent  term  is  small.  If  a  substance  could  be  synthesized  or 
designed  such  that 
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then  there  woulu  be  no  aberration  for  that  polarization  in  the  quartic  approximation 
to  the  truncated  Gaussian  beam.  In  fact,  this  effect  is  quite  prominent  in  the 
I  -  VII  and  II  -  VII  compounds  shown  in  Tables  1  and  2.  This  shows  that  it  may  be 
possible  to  find  or  design  An  ionic  material  where  Eq.  (57)  is  fulfilled.  Satisfying 
Eq.  (57)  would  not  eliminate  all  aberration,  since  there  would  still  be  a  small 
amount  resulting  from  the  terms  that  we  neglected  in  Eq.  (31). 

6.  1.  3  COMPOSITE  WINDOW 

It  is  possible,  in  principle,  to  find  two  different  materials,  one  with  a  positive 
0w  and  the  other  with  a  negative  8W  of  the  same  magnitude,  such  that  when  they 
act  in  series  as  a  composite  window  the  aberrations  nearly  cancel  out.  If  two 
such  materials  could  be  found,  a  diffraction-limited  focus  could  be  maintained  for 
a  long  time. 

6. 2  Engineering-Type  Solutions 

6.2.1  COOLING 

One  method  of  minimizing  the  thermally-induced  aberrations  is  to  cool  the 
window  in  such  a  way  as  to  minimize  the  temperature  gradient  in  the  p  direction. 
We  can  imagine  several  schemes  of  forced  face-cooling  below  ambient  tempera¬ 
ture  which  may  accomplish  this,  and  these,  as  well  as  other  cooling  configurations, 
should  be  thoroughly  studied.  If  one  face-cools  too  rapidly,  however,  large 
temperature  gradients  in  the  z  direction  will  be  set  up  which  could  either  per¬ 
manently  strain  the  window  or  crack  it. 

6.2.2  MOVING  WINDOW  OR  MOVING  BEAM 

Another  engineering  possibility  suggested  by  M.  Sparks  is  to  minimize  ther¬ 
mal  effects  by  moving  a  large  window  with  respect  to  a  smaller  laser  beam.  This 
would  allow  the  beam  to  travel  through  fresher  material  than  if  the  entire  system 
had  been  fixed.  After  the  exposed  portion  of  the  window  rotated  away  from  the 
beam,  it  could  be  cooled. 

6. 2.  3  OPTICAL  COMPENSATION 

A  very  appealing  solution  would  ba  to  compensate  fully  for  window  aberrations 
by  designing  a  continually  changing  fore-optical  system.  This  system  would  have 
to  change  shape  in  such  a  manner  that  it  caused  the  incident  angle  at  the  window  to 
be  given  by: 
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When  9  comp  was  added  to  9i  I  Siven  by  Eq.  (49)] ,  the  resulting  total  angle  of 
deviation  would  be: 

®tot  =  -lf*  (59) 

o 

This  means  that  no  spherical  aberration  would  occur  and  that  the  beam  would 
remain  diffraction-limited  for  a  long  time. 

7.  PROGRAM  RESEARCH  AND  DEVELOPMENT  NEEDS 

In  this  section,  we  list  some  laser  window  program  needs  in  the  area  of 
materials  research  and  development.  *  We  divide  them  into  three  areas:  im¬ 
mediate  R  and  D  objectives  as  well  as  longer-term  experimental  and  theoretical 
studies.  ** 

7.1  Some  Specific  Immediate  R  and  D  Objectives 
7.1.1  TO  OBTAIN  DATA 

More  and  better  data  on  key  material  parameters  are  needed  at  both  10.  6/u 
and  3  to  5u  so  that  prospective  window  candidates  may  be  properly  evaluated. 
These  important  parameters  are: 

(a)  bulk  absorption  coefficient  0A  ; 

(b)  specific  heat  c'  ; 

(c)  linear  expansion  coefficient  a  : 

(d)  index  of  refraction  n  : 

(e)  temperature  coefficient  of  index  of  refraction  ( —p  ^ 

WT/c=o 

(f )  Poisson's  ratio  v  ; 

(g)  Pockel's  elasto-optic  constants  p„;  and, 

(h)  thermal  conductivity  K  . 

*  A  brief  discussion  of  some  of  these  needs  was  presented  by  Dr.  J.  R. 
Jasperse  at  the  ARPA-sponsored  meeting  on  laser  window  materials  at  Wood's 
Hole,  Mass,  in  July,  1971. 

**  This  list  should  not  be  considered  as  exhaustive  or  as  representative  of 
the  current  laser  window  program  at  AFCRL,  It  gives  the  authors'  thoughts  on 
some  of  the  more  important  program  needs  as  of  July  1971. 
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For  many  materials  of  interest,  data  are  just  not  available;  whereas  in  other 
cases  they  are  not  reliable. 

7. 1. 2  TO  IMPROVE  AND  DEVELOP  SOME  SPECIFIC  CANDIDATE 

MATERIALS 

(a)  Alkali  Halides 

Here,  work  in  the  immediate  future  should  be  centered  around 
the  two  objectives  of  improving  the  mechanical  properties  of  the  I  -  VII's 
and  reducing  their  propensity  to  pick  up  water. 

(b)  II  -  Vi's  and  III  -  V's 

Two  aspects  of  research  and  development  should  be  emphasized 
here:  (1)  growth  of  crystals  of  increased  purity  and  improved  optical 
quality,  so  as  to  produce  materials  that  achieve  a  0^  approaching  its 
intrinsic  value;  and,  (2)  development  of  large-size  growth  techniques 
since  difficulties  are  anticipated  in  growing  large  crystals  of  these 
materials. 

(c)  Glasses 

Research  and  development  on  glasses  should  be  geared  toward 
solving  the  problem  of  designing  a  glass  with  a  small  0A  and  with  an 
increased  thermal  conductivity.  Large-size  window  fabrication  should 
not  be  a  major  problem  for  the  glasses,  since  modest-sized  castings 
are  currently  available  with  present  technology. 

(d)  Polycrystalline  Materials 

The  important  objective  here  is  to  make  polycrystalline  aggregates 
with  0^  close  to  intrinsic  values.  This  would  provide  window  materials 
which  could  be  formed  in  large  sizes  and  yet  have  small  values  for  0^  . 

7.2  Some  Specific  Experimental  Studies  Needed 

(1)  The  optical  properties  of  polycrystalline  materials  as  a  function  of  grain 
size.  These  studies  should  be  carried  out  for  both  a  cubic  and  a  noncubic  ma¬ 
terial  chosen  such  that  grain  sizes  can  be  varied  over  a  wide  range  compared 
with  the  wavelength  of  the  light. 

(2)  The  influence  of  impurities  and  defect  structure  on  the  optical  properties 
of  candidate  crystals. 

(3)  The  surface  and  bulk  scattering  properties  of  real  candidate  materials. 

(4)  The  surface  and  bulk  absorption  properties  of  real  candidate  materials. 

(5)  Determination  of  the  intrinsic  dielectric  breakdown  limit  in  real  candidate 
materials  as  power  densities  are  increased  over  those  currently  anticipated. 
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(6)  The  development  of  coating  technology  for  all  candidate  materials, 
especially  for  the  alkali  halides  where  difficulties  are  known  to  exist. 

(7)  Energy-loss  mechanisms  studies. 

7.3  Some  Specific  Theoretical  Studies  Needed 

(1)  A  fundamental  theory  for  n  and  (3^  in  dielectric  crystals  which  is  accurate 
in  the  vicinity  of  10.  6  microns. 

(2)  A  fundamental  study  of  why  (■§■%)  is  negative  for  ionic  crystals  and 

' o-o 

positive  for  covalent  crystals. 

(3)  A  fundamental  model  for  the  elasto-optic  coefficients  accurate  in  the 
vicinity  of  10.  6  microns. 

(4)  A  study  of  the  effect  of  coherent  light  on  the  absorption  of  power  in  a 
disc-shaped  laser  window. 

(5)  A  study  of  the  effect  of  surface  absorption  on  thermal  lensing  and  on 
localized  heating. 

(6)  A  study  of  the  effect  of  surface  and  bulk  scattering  on  the  temperature 
distribution  in  the  laser  window. 

(7)  Studies  on  how  to  design  a  new  material  with  a  maximum  diffraction- 
limited  time  and  witn  other  desirable  physical  properties. 

(8)  Calculations  of  u  versus  k  in  prospective  window  materials  to  include 
some  ternary  compounds. 
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Appendix  A 


Derivation  of  At  Polorisation-Dtptndtnt  Refractive  Indicts 

as  Powor  Series  in  p/p0 


In  order  to  evaluate  the  two  integrals  in  Eqs.  (19a  and  b),  we  must  first 
determine  AT.  Substituting  Eqs.  (2)  and  (3)  into  (16)  and  replacing  p/pQ  by  q,  we 
get: 


AT(q,t)  =  g(t)  2  T_.  q^ 


Then, 


fq dx  x  AT(x.t)  =  g(t)  ^2  T2<  p  dx  x2/+1  =g(t)  ^2  -^f§q2l+2 


/~o  2<+2 


f1  dxx  ATfe.t)  =  g(t)  2 
t  =c 


«  T„ 


When  these  integrals  are  inserted  into  Eqs.  (19),  there  results: 


a  =«Y 
PP 


(A-4) 


{[pu-21'I>12][1-‘'2<]+  [O-IPk-Ph]  [l-(21+l)q2(]  | , 
This  can  be  rearranged  to: 


(A-7) 


np  =  n  +  g(t) 


jg ?0T2^  | ( 9T  )a=Q  q2<  "  2(21+2)  [(1_,/)pll 

+  (l-3v)p12  -q2<  |(l-(2f+l)v)  pn  +  (2l+l-(2*+3)i/)p12jj  J  . 


(A-8) 


We  now  partition  the  summation  into  two  terms:  one  being  the  summand  evaluated 
at  I  =  o,  the  other  being  the  sum  from  one  to  infinity.  When  this  is  done,  one  can 
easily  see  that  there  will  be  some  terms  that  do  not  contain  q.  These  are  collected 


under  the  symbol  np(0,  t),  leaving: 


.np(0..)  +  g(0(E  t2j  q2' {(I%)ct_o 


jSXL 


2(21+2) 


[(l-(2<+l)v)p11+  (21+1  -  (2f+3)v)p12]|. 


(A-9) 
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This  can  be  rewritten  as: 

np=np<°'t)  +  g(t)i?i  n2i  S'  ‘  (A- 10) 

which  is  Eq.  (20aland  where 

"2<  *  T2I  [(&r)ff.0  +  +  (2<tl-(2/+3Wp12]  |  .  (A-ID 

which  is  Eq.  (21a). 

The  same  line  of  reasoning  pertains  to  n^(q,  t).  In  its  counterpart  to  Eq.  (A-6), 
the  np  is  replaced  by  n^.  and  the  CQ  and  Cj  are  interchanged.  Thus,  the  counter¬ 
part  to  what  is  inside  the  braced  term  in  Eq.  (A-7)  becomes: 

| [a-vfeja ^2t]  +  [pii - 2vPn] [  1  -  <2<+1> q2< ]  } * 

which  can  be  rearranged  to: 

[(l-vfcjj  +  (l-3v)p12  -  q2<  ^f+l-l/Jpjj  +  (l-(4f +3)v)p12^J  . 

\gain,  partitioning  the  sum  into  two  terms  and  collecting  all  terms  not  containing 
q  under  the  symbol  n^(0,  t),  this  leaves: 

n*  =  V0*  *>  +  S  » J4  q21  .  (A- 12) 

which  is  Eq.  (20b),  where 

°2I  *  T2l|(lr)a=0  +  2§m)  [»l+l-»h»n  +  (l-(4/+3)v)p12J  \  .  (A- 13) 


which  is  Eq.  (2  lb) . 


Appendix  B 

Approximate  Solution  to  Roy  Tract  Equation 


According  to  the  ray-trace  formulas,  the  radius  of  curvature  R  at  any  point 

c 

on  the  ray  trajectory  is  (Born  et  al,  1964): 

*  P.  V  (In  n)  *  |v(f  n  n)|  cos  0  ,  (B-l) 

c 

where  &  is  the  unit  normal  along  Rfi  and  n  represents  either  np  or  n^  (See  Figure  1). 
From  elementary  calculus  the  radius  of  curvature  of  a  simple  plane  curve  is 


Consider  an  index  of  refraction  which  depends  only  on  p,  and  which  is  either 
monotonically  increasing  or  monotonically  decreasing,  so  that  the  algebraic  sign 
of  Rc  does  not  change.  If  we  equate  the  right-hand  side  of  the  two  equations  above, 
and  replace  cos  6  by  ( 1  +  tan2®]  ~^2,  where  tan  9  =  dp/dz,  then  we  obtain  the 
following  nonlinear,  second-order  differential  equation  for  the  ray  trajectory: 


Preceding  page  blank 


Sw 


’■v"  *  v»* ^ tJ531T^.  ¥  ■‘V".  .¥  i* - i*jH*7S5>V»- f t *  ^j; 


?  mJSWJwnw)^^ 


For  small  angles  of  deviation,  ('"j^)2  is  negligible  compared  to  unity.  * 

Since we  are  assuming  refractive  indices  of  the  form  shown  in  Eqs.  (20a)  or 
(20b),  the  dn/dp  are  calculable  as  power  series  and  the  n  in  the  denominator  of 
Eq.  (B-2)  can  be  approximated  by  n(0,t).  The  nonlinear  equation  becomes: 


-4-  (^-)  ?  "21  *  01“/"' 

•3z2  yi(0,t>poy  1*1  21  ° 


Looking  for  a  solution  of  the  form 


p(z)  *  pj  (1  +  y(z)  J  . 
where  |y|  «  1  ,  we  obtain 

“ M*(^j  *  y>2"' • « ■ 

21-1 

Expanding  (1  +  y)  up  to  first-order  in  y,  the  differential  equation  can  be 
written  as: 


2  K(t)  -  2  K'(t)  y  =  0  . 
dz£ 


where 


K,t)-  MCMJp ,J,  "2l'(p‘  ) 


K'ft>  •  n(0.(iipop.  "21  '  «' -1>(p,/p0>"'1  • 

Employing  the  boundary  conditions  that  y  =  dy/dz  =  0  at  z  =  0,  we  can  get 
2  possible  solutions  for  y,  depending  on  the  sign  of  K'.  For  K'  >  0, 


♦This  is  equivalent  to  the  paraxial-ray  approximation. 
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y  *-||t  |coah»/2K'  z  -  1^. 

3  6 

Or,  expanding  up  to  order  K  z  ,  remembering  that  K  and  K1  have  the  same  order 
of  magnitude: 

•  ySKx2+-y!  z4  +  o{k3z6 
whereas,  for  K'  <  0  , 


Again,  expanding. 


Thus,  for  either  case,  the  solution  can  be  expressed  as: 

p(z)  -  Pt  [l  +  Kz2  +  bJ-Ki  z4  +  o{  |K|3  z6  }].  (B-8) 

where  K  can  be  either  positive  or  negative. 


Appendix  D 

The  Circle  of  Leost  Confusion  From  Goomotric  Optics 


In  this  Appendix  we  calculate  the  radius  of  the  circle  of  least  confusion  and  its 
distance  from  the  window  using  geometric  optics  for  either  polarization.  We 
define  the  distance  from  the  window  to  the  point  where  the  ray  crosses  the  z  axis 
by: 


where. 


7<p/p0>=  -®tot' 

in  the  small  angle-of-deviation  approximation.  Employing  Eq.  (51)  in  the  above: 


*■(£)[= -f*(£)2} 


-"■{**£  i* -p.  ft/]  r 


Since  RQ  is  always  in  front  of  the  window  (i.  e. ,  positive),  we  will  confine  our 
attention  only  to  those  circumstances  in  which  R(p)  also  remains  positive. 
Regardless  of  the  magnitudes  and  signs  of  Fj  and  F2,  this  condition  can  always 
be  satisfied  provided  1  -  Fg  (p/pQ)2  is  always  greater  than  -  Pq/Fj  Rq  •  Then  R, 
and  henct  0^,  rema*n  monotonic  with  p  . 

In  Figure  5,  we  depict  two  of  the  rays  for  the  case  in  which  R  increases 
monotonically  with  p.  The  ray  emanating  from  p  =  pQ  makes  an  angle  of  yQ  with 
the  horizontal  and  has  an  equation  given  by: 

p =  ~  y oz  +  po  *  *D"5 

A  second  ray,  emanating  from  p  =  -p',  is  shown  coming  from  the  bottom  half  of 
the  window  for  reasons  that  will  soon  become  obvious.  It  makes  an  angle  of  7' 
with  the  horizontal  so  that  its  equation  is: 


p  =  7'  z  -  p' 


The  two  rays  intersect  at  a  distance  r  from  the  origin  given  by: 


="r^7r  [Po7'  ’  P  7o]* 


The  maz.imum  value  of  r,  which  we  call  rc>  is  the  radius  of  the  circle  of  least 
confusion.  To  obtain  rc,  we  first  maximize  r  in  Eq.  (D-7)  with  respect  to  p'/p0 
and  get  a  corresponding  value  for  p'.  This  particular  p\  together  with  pQ,  are  all 
that  is  necessary  to  prescribe  those  two  rays  which  are  sufficient  to  determine  the 
circle  of  least  confusion. 

Taking  the  derivative  of  Eq.  (D-7)  with  respect  to  p'/p0  results  in: 

'  *  h  (y  y°)  ■  ®'8) 


Setting  this  equal  to  zero  yields: 


(iL  +  l\  -  d^'  =  7'  + 

\po  ;d(P'/pj 


Now  utilizing  Eq.  (D-3)  at  both  pQ  and  p',  the  above  reduces  to: 


2(£) + ’(-fe)2  * 1  ■ 0  • 


(D- 10) 
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Of  the  three  solutions,  namely,  1/2,  -1,  and  -1,  only  p'/p0  =  1/2  is  the  phys:o*»lly 
acceptable  solution.  Substituting  this  into  (D-7)  and  changing  r  to  rc,  the  radius 
of  the  circle  of  least  confusion,  for  either  polarization,  is  now  given  as: 


(D-ll) 


The  distance  of  the  circle  of  least  confusion  from  the  window,  can  be 
obtained  from  either  Eqs.  (D-5)  or  (D-6)  if  z  is  replaced  by  d£  and  p  by  r^  . 
Solving  for  results  in: 


When  R(p)  is  a  monotonically  decreasing  function  of  p,  we  find  that  the  same 

equations  as  given  above  hold  and  that  r  and  have  the  same  forms  as  in 

c  c 

Eqs.  (D-ll  and  12).  The  signs  of  the  F-  functions  may  differ,  however,  and  only 
the  magnitude  of  r  j,  is  of  consequence. 


